Astronomy <fc Astrophysics manuscript no. 1219 


February 2, 2008 


(DOI: will be inserted by hand later) 





Statistical properties of exoplanets IV. The period-eccentricity 
relations of exoplanets and of binary stars 

J.L. Halbwachs 1 , M. Mayor 2 , and S. Udry 2 

1 Observatoire Astronomique de Strasbourg (UMR 7550), 11 rue de l'Universite, F-67 000 Strasbourg, France 
email: halbwachs@astro.u-strasbg.fr 

2 Geneva Observatory, 51 chemin des Maillettes, CH-1290 Sauverny, Switzerland 
email: Michel. Mayor@obs.unige.ch, Stephane.Udry@obs.unige.ch 

Received 3 May 2004; accepted 19 October 2004 

Abstract. A sample of spectroscopic binaries and a sample of single planetary systems, both having main-sequence 
solar-type primary components, are selected in order to compare their eccentricities. The positions of the objects 
in the (P-(l — e 2 ) 3//2 , e) plane is used to determine parts in the period-eccentricity diagram that are not affected 
by tidal circularization. The original eccentricities of binaries and planets are derived and compared. They seem 
to be weakly or not at all correlated with period in both samples, but two major differences are found : 

(1) The tidal circularization of planetary orbits is almost complete for periods shorter than 5 days, but it is not 
visible when P.(l — e 2 ) 3//2 is longer than this limit. This suggests that the circularization occurs rapidly after 
the end of the migration process and is probably simultaneous with the end of the formation of the planet. By 
contrast, we confirm that the circularization of the binary orbits is a process still progressing a long time after 
the formation of the systems. 

(2) Beyond the circularization limit, the eccentricities of the orbits of the planets are significantly smaller than 
those of binary orbits, and this discrepancy cannot be due to a selection effect. Moreover, the eccentricities of 
binaries with small mass ratios are quite similar to those of all binaries with q < 0.8. This suggests that the low 
eccentricities of exoplanet orbits are not a consequence of low-mass secondaries in a universal process. 

These remarks are in favor of the idea that binaries and exoplanets are two different classes of object from the 
point of view of their formation. 
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systems: formation 



1. Introduction 

It is well known that the orbits of the exoplanets with 
periods larger than 5 or 6 days have eccentricities signif- 
icantly larger than those of the giant planets of the solar 
system. Several mechanisms were proposed to explain this 
feature, but, up to now, none is fully convincing. It was 
proposed that eccentric orbits could be a consequence of 
the dynamic evolution of systems initially involving sev- 
eral planets (Rasio & Ford lTOOol Lin & Ida lTOOTl Ford et 
al. EMI Papaloizou & Terauem I200T1 Rice et al. 12005)) . 
but these models fail to produce the frequency of giant 
planets with semi-major axes smaller than about 1 AU. 
The giant planets close to their harboring stars are of- 
ten assumed to be produced by migration within a disk 
(Ward 119971 Masset & Papaloizou 2003) and references 
therein), but this process hardly produces eccentric or- 
bits (Papaloizou et al. 120011 Thommes & Lissauer l2003[l . 
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although Goldreich & Sari l)200"5)l and Woolfson (2003) 
leave some room for hope. Therefore, it is tempting to 
consider that the exoplanets are generated by the same 
process as binary stars (Stepinski & Black 2000). This 
implies that giant exoplanets were not formed by gas ac- 
cretion onto a heavy rocky core, as usually assumed, but 
by an alternative process. They could come from disk in- 
stabilities (Mayer et al. 120021 Boss l2002ll2005)l . but inward 
migration in a disk is then invoked again to explain the 
short-period orbits; alternatively, planets could be gener- 
ated by fragmentation of a collapsing protostellar cloud, 
via filament condensation and capture (Oxley & Woolfson 
2004}, or even exactly as stellar components in binary sys- 
tems (see discussion in Bodenheimer et al. 120001 and refer- 
ences therein). However, these models may be efficient in 
forming massive planets or brown dwarfs, but not planets 
around 1 Jupiter mass or less. 

Note that the binary formation models are not very 
satisfactory either (see the review by Tohline l2()()2|l . The 
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large eccentricities of binaries are explained by fragmen- 
tation of collapsing cores and subsequent interactions be- 
tween the forming stars (Bate et al. 120021 Goodwin et al. 
2004), but, as for exoplanets, the simulations do not pro- 
vide the high frequency of close systems. Moreover, statis- 
tical investigations on main sequence binaries (Halbwachs 
et al. 120031 Paper I hereafter) have shown that the close bi- 
naries (i.e. with semi-major axes less than a few AU) con- 
sist in two populations : one with large eccentricities and 
mass ratios less than 0.8 ("non-twins" hereafter), and one 
with moderate eccentricities and nearly identical compo- 
nents ( "twins" ) . Additionally, the twins are more frequent 
among short-period binaries than among the others. At 
first, these properties were derived from binaries with F7- 
K primary components, but they are also valid for M-type 
dwarfs (Marchal et al. l2UU3|l . 

In the present paper, the period-eccentricity diagram 
is used to compare the exoplanets with the binary stars : 
our main purpose is to investigate if the properties of exo- 
planets may be considered as an extrapolation of the prop- 
erties of binaries in the range of very low mass ratios. This 
would indicate whether the formation processes of these 
objects are similar. In the course of the paper, a few other 
points are also treated : (1) the correlation between the 
eccentricity and the period or the angular momentum, (2) 
the relation between the eccentricity and the metallicity of 
planets, (3) the original distributions of eccentricities for 
planets and for binaries, considering the twins separately. 
Comparisons between planets and binaries in the period- 
eccentricity diagram were already presented by Mayor et 
al. (EOOIl), Mazeh & Zucker pTtTT)) . and Udry (JSHU), who 
concluded that planets and binaries are very similar when 
periods longer than 50 days are considered. However, their 
samples contained around 30 or 40 planets, and a many 
others have been discovered since. The question needs to 
be re-considered. 

The interpretation of the period-eccentricity diagram 
is rather complex, and our investigations are based on the 
method presented in Sect. [5] Sect. 01 is dedicated to the 
binaries; we investigate if, additionally to the twins, other 
classes of mass ratio have specific distributions of eccen- 
tricities. A similar treatment is applied to exoplanets in 
Sect. 0| Binaries and exoplanets are compared in Sect. [5] 
in which we pay attention to the difference in the selection 
effects of both samples. The consequences of our results 
are discussed in Sect. El 

2. Method 

2.1. Tidal effects 

We must pay attention to the fact that the periods (P) and 
the eccentricities (e) of the objects are modified by tidal 
interactions, especially when P is short. As a consequence, 
the (P - e) diagram may schematically be divided into 
two parts : the short periods, where the orbits are circu- 
lar or have low eccentricities, and the periods longer than 
the circularization limit, hereafter called P cu toff (Mayor & 



Mermilliod lTggl Duquennoy & Mavor lTMTjl . Several the- 
oretical models were proposed to derived P cu toff) and the 
treatment is not the same for binaries (Zahn 1992, Hut 
1981, 1982, Keppens ll997l and references therein) and for 
planets (Goldreich & Soter HMHl TrillingEIIIH) ■ Moreover, 
several physical processes are invoked, each of them de- 
pending differently on the mass ratios of the systems. 

Despite the complexity of the process, a few simple 
guidelines may be drawn. First of all, the efficiency of tides 
in modifying the orbits is very sensitive to the distance be- 
tween the components. For a given system, the tidal torque 
depends on the orientation of the tidal bulge and on the 
separation between the components, r. It varies as 1/r 6 
(Lecar et al. I1976JI . Therefore, for systems differing only 
by period, the transition from circularized orbits to orbits 
practically unaffected by tides is a narrow strip in the (P 
- e) diagram (see the simulations by Witte & Savonije, 
2002). However, this does not mean that the systems with 
P > Pcutoff may have any eccentricity. For a given period 
the systems with eccentric orbits have components much 
closer than the semi-major axis during a part of the pe- 
riod, and below a certain limit the orbit rapidly becomes 
circular. As a consequence, the upper part of the (P - e) 
diagram is cleared even for P > P cu toff- Note that the or- 
bits that were originally eccentric do not keep the same 
period when they are evolving towards e = : when the 
primary star is a slow rotator, the orbit is circularized 
keeping the orbital angular momentum unchanged (Witte 
& SavonijeHinS Hurley et al. lMJ2) . Therefore, the semi- 
latus rectum, r sr = a(l — e 2 ), is conserved and it becomes 
the radius of the final circular orbit. As a consequence, 
the final period is : 

P sr = P(l- e 2 ) 3 / 2 (1) 

where P and e refer to the original state of the sys- 
tem. This gives us a simple but efficient way to explore 
the transition from circularized orbits to orbits unaffected 
by tidal effects. In place of a (P - e) diagram, the ob- 
jects are plotted in the (P sr - e) plane. Therefore, the 
evolution path toward a circular orbit is a vertical line 
in the diagram. Moreover, note that, for a wide range of 
eccentricities, the mean value of 1/r 6 for a complete or- 
bit is approximately l/r sr 6 (Fig. QJ. As a consequence, 
we expect that, when r sr is small enough to permit effi- 
cient tidal effects for a given eccentricity, these effects will 
remain important during the evolution of the orbit, until 
they eventually lead to circularization. Therefore, the bor- 
der between the circularized systems and the area nearly 
unaffected by tides in the (P sr - e) diagram should appear 
very clearly. If the circularization were not at all related 
to the ages of the systems, the diagram should show a 
strong contrast, with the circular orbits on the left-hand 
side, and the orbits with any eccentricities, including the 
largest ones, immediately beyond the circularization limit. 
This is almost what is observed in reality, especially for 
the exoplanets (see Fig. [3] and 0] below). 

In practice, however, a real sample is selected up to a 
maximum period Pviax, and not all the eccentricities are 
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Fig. 1. The mean value of the 1/r 6 term of the tidal 
torque, compared to l/r„ , where r sr is the semilatus 
rectum. Except for very large eccentricities, l/r e T is an 
acceptable approximation of < 1/r 6 >. 

permitted in the (P sr - e) diagram. For a given P sr , the 
eccentricities range from to the limit : 

e S up = ^l-(PsrAPMax) 2 / 3 . (2) 

With Pcutoff so determined, we derive from Eq. (JJJ 
the maximum eccentricity that the systems unaffected by 
tides in the (P - e) diagram may have : 

e Max = yjl- (Pcutoff /P) 2 / 3 - (3) 

For comparing two samples which were have been dif- 
ferently affected by tidal circularization, it is necessary to 
restrict the comparison to the smallest of the two limits in 
eMax in order to find discrepancies coming only from the 
original distributions of e. 

2.2. Comparison of two samples 

The comparison between different samples of orbiting 
systems is based on the median eccentricities. Two ap- 
proaches are used : the first is a visual examination, and 
the second is a statistical test. In the first approach, the 
median eccentricity of each sample is derived as a func- 
tion of the period. In practice, for any period Pi > P cu toff, 
the median eccentricity is derived from 6 systems taken 
among the periods closest to Pj. Except near the limit 
of the period range, 3 of these 6 systems have P < Pi 
and the 3 others have P > Pj. The second approach is 
the statistical test used in Paper I : the two samples are 
merged, and the range of periods longer than P cu toff is di- 
vided into bins, each containing 12 systems, except for the 
last one which may contain up to 23 systems. The com- 
mon median eccentricity is derived in each bin, and the 
systems below the median are counted for one of the two 
samples. If this sample actually belongs to the same sta- 
tistical population as the other, the probability of getting 
any count, P(k), obeys the hypergeometric distribution. 



The rejection threshold of the null hypothesis, Hq : "all 
the systems are equivalent from the point of view of the 
eccentricities" is then derived. When the count k is less 
than half the population of the considered sample, the re- 
jection threshold of H in a two-sided test is 2 x P(i < k); 
on the contrary, it is 2 x P(i > k) when k is larger than 
the expected number. 

Note that setting the content of the bins to 12 systems 
is a bit arbitrary. This number is neither too small to 
derive a reliable median nor too large to have a nearly 
constant period distribution within each bin. Using other 
numbers close to 12 would give other results, but it was 
verified that the differences are not important. 

3. The binaries 

3.1. The binary sample 

The so-called extended sample of F7-K dwarf binaries se- 
lected in Paper I is used again. It consists of 89 spectro- 
scopic binaries (SB) found in the solar neighborhood or 
in open clusters, with periods of up to 10 years. We al- 
ready know that the twins have, on average, eccentricities 
smaller than the other binaries. However, before compar- 
ing binaries to exoplanets, it is worthwhile to sec if the 
eccentricities of non-twin binaries depend on the mass ra- 
tios. 

The mass ratios q = M.2I M.\ of the SB in the sample 
have been fixed for 58 binaries, thanks to the combina- 
tion of the SB orbital elements with Hipparcos astrometric 
observations, or with photometric sequences in the open 
clusters (Paper I). For the other 31 SB, we derive intervals 
containing the actual mass ratios. The minimum limits 
are computed from the mass functions; the maxima are 
obtained differently for the nearby SB and for the cluster 
SB : the former all have q < 0.65, since otherwise they 
would be double-lined SB with known q, and the latter 
have limits coming from their positions in the photometric 
sequence of the cluster. 

In order to make visible a possible relation between the 
mass ratios and the eccentricities, the SB are distributed 
in several groups : q < 0.40 (16 SB), 0.40 < q < 0.80 (20 
SB), and twins (27 SB); we still add a group containing 
all the SB with q < 0.80 (62 SB, including the 36 already 
in the first two groups). 

3.2. Limit of tidal circularization 

The SB are plotted in the (P sr - e) diagram (Fig. 0), to 
delimit the part of the diagram affected by tidal circular- 
ization. 

Several relevant features appear in Fig. [21 First of all, 
the range of eccentricities jumps from to almost 1 be- 
tween 4 and 9 days. All SB with P sr shorter than 4 days 
are at present on circular orbits, and no evidence of cir- 
cularization is visible when P sr exceeds 10 days. The SB 
with P sr between 4 and around 10 days have medium ec- 
centricities, and some of them are even circularized. 
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Fig. 2. Distribution of the SB in the (P sr - e) diagram, 
where P sr is the period that the SB would have if their 
orbits became circular with the same angular momenta. 
The symbols represent the systems with different mass 
ratios. The code is as follows : q < 0.4, circles; 0.4 < q < 
0.8, squares; q < 0.8 (maximum q for SB1) and minimum 
q < 0.4, small diamonds; q > 0.8 (twins), open triangles. 
The thin dotted line is the limit es U p corresponding to a 
maximum period of 10 years, as derived from Eq. PJl. 




Fig. 3. Distribution of the SB in the period-eccentricity 
diagram. The symbols are the same as in Fig. [21 The lines 
refer to the median eccentricities of the different classes 
of mass ratios : q < 0.4, thick full line; 0.4 < q < 0.8, 
thick dashes; all q < 0.8, thick dot-dashed line; q > 0.8 
(twins), thin line. The 2 thin dotted lines are the limits 
of the location of SB with P sr between 5 and 10 days, as 
derived from Eq. ©. 



At least one of the short period SB may have been 
generated with a very large eccentricity, and may be too 
young for to have been circularized : KW 181 has the 
most eccentric orbit among the periods shorter than 10 
days. This system belongs to the Praesepe cluster, and its 
age is therefore only 0.8 10 9 years. Its circularization will 
be complete within only 0.1 10 9 years (Duquennoy et al. 

The other two SB with e > 0.1 and P > 10 days 
may have eccentricities due to perturbations by a third 
component (Kozai mechanism or secular perturbations, 
see Mazeh & Shaham HHTHI : GJ 719 has a CPM compan- 
ion with a minimum separation of 280 AU (Zuckerman et 
al. ©97), and GJ 233 is a visual binary with a possible 
period of 200 years (Heintz Il988[) . When these 3 SB are 
discarded, we find P cu toff ~ 8 or 10 days, as in Duquennoy 
& Mayor HTWd . 

Nevertheless, it is striking that all the SB with ec- 
centric orbits and P < 10 days have large mass ratios. 
In order to see if this feature is significant we look at 
the period-eccentricity diagram of the 205 SB found by 
Latham et al. I|2002|l among stars with large proper mo- 
tions. They also found a range of period where circular 
orbits and eccentric orbits both exist, but between around 
10 and 20 days (see their Fig. 9). It is visible in their plot 
that the double- lined SB (SB2), which have the largest 
q, are not abnormally frequent among the systems with 
eccentric orbits and short periods. Therefore, we admit 
that the frequency of twins with large eccentricity and 
P < 10 days is just due to chance. 

The sample of Latham et al. contains a few SB with cir- 
cular orbits and periods between 10 and 20 days. However, 
they are stars with large proper motions, and they gener- 



ally belong to the old galactic disk or even to the halo. 
Therefore, the long periods of some circularized orbits 
may be an effect of the ages of the systems, in agree- 
ment with theoretical predictions (see e.g. Duquennoy et 
al. H992f) . Moreover, these old stars are not representative 
of the stars observed for planet detection. For the par- 
ent population of stars harboring planets, the limit of the 
area affected by tidal effects in the (P SI - e) diagram is 
Pcutoff = 10 days. 

3.3. The binaries in the (P - e) diagram 

The (P - e) diagram of the SB is plotted in Fig. The 
values of the median of the four classes of q are drawn in 
this figure for visual comparison. 

As shown in Paper I, the twins often have below- 
average eccentricities when periods longer than 5 days are 
considered; the probability to get so large a discrepancy 
by chance, as derived by the two-sided test, is only 2.7 %. 
By contrast, the two groups with q < 0.80 have nearly the 
same distribution of eccentricities : the significance of the 
two-sided test is 100 %. In order to check that even the SB 
with the smallest mass ratios have the same eccentricity 
distribution as the others, the test is done again by com- 
paring the SB with q < 0.25 (10 SB with P > 10 days) 
to those with 0.25 < q < 0.80 (32 SB with P > 10 days). 
Again, exactly half of the low-mass ratio SB are below the 
median, providing a 100 % significance. Therefore, all the 
non-twin SB may be considered together in the compari- 
son to the exoplanets. 

We now come back to the distribution of e of all SB. It 
is clearly visible in Fig.|21that the median eccentricity in- 
creases with the period. However, we suspect that this may 
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be entirely explained by tidal circularization, because the 
orbits having initially P > P cu toff and e > eMax are now 
circular with P = P sr < P cu toff- I n other words, we want 
to see if the distribution of eccentricities depends on the 
periods, apart from the cut-off at e = eMax- A Spearman 
test is performed to check this hypothesis. In order to dis- 
card the area affected by tidal circularization, the test is 
restricted to the rectangular box (P > 20 d, e < 0.61). 
The non-twins and the twins are considered separately. 
The Spearman correlation coefficient is 0.26 for the for- 
mer, and 0.29 for the latter; taking into account the num- 
bers of objects, the probabilities to get by chance values 
so far from zero are 10 % and 20 % respectively. Although 
these levels of significance are a bit low, they are still too 
large to reject the hypothesis that the eccentricities are 
not correlated with the periods as soon as the semilatus 
rectum of the orbit is larger than the radius of a circular 
orbit with period P cu toff- 

3.4. Correlation eccentricity versus angular momentum 

Since P sr is related to the angular momentum of the orbit, 
it seems relevant to see if the eccentricities are correlated 
with this parameter. This question looks similar to the (P 

- e) correlation investigated just above, but it is different 
in reality, since changing the (P - e) plane into the (P sr - 
e) plane modifies the density in relation with the distribu- 
tion of the periods. Therefore, the absence of correlation 
between P and e does not necessarily imply an absence of 
correlation between P sr and e, and vice versa. 

A Spearman test is used again, but in a box in the (P sr 

- e) diagram. Since the selection of the sample was limited 
by the condition P < 10 years, we must take into account 
the limit es up derived from Eq. (0) (the thin dotted line 
in Fig. 0). Therefore, the limits of the box considered in 
the Spearman test are P sr £ [10 d, 789 d] and e < 0.8, in 
order to avoid the area with e > es U p- Again, the rejection 
threshold of the hypothesis that e and P sr are correlated 
is between 10 and 20 %, and the eccentricities may not 
depend on the angular momenta of the SB. 

4. The exoplanets 

4.1. The planet sample 

We start from the up-to-date list of exoplanets that is 
provided on the Geneva web site 1 . However, several plan- 
ets cannot be used in a comparison to the SB, for several 
reasons : 

— Since the binaries all have main-sequence primary 
components, the planets orbiting subgiant or giant 
stars are discarded. The specific problem of tidal cir- 
cularization in a system containing an evolved star is 
thus avoided, as well as the uncertainty coming from 
the evaluation of the mass of an evolved primary com- 
ponent. 

http:/ /obswww.unige.ch/Exoplanets 



— The eccentricities must be reliable, and therefore de- 
rived from an orbit with good quality. For that pur- 
pose, the residuals of the radial velocity (RV) measure- 
ments, rms, are compared to the semi-amplitudes of 
the spectroscopic orbit, K. Five planets having rms 
larger than if/3 are discarded. This criterion looks 
a bit rough, since several other factors could also be 
taken into account, such as the number of measure- 
ments, the phase distribution, and the fact that some 
orbits refer to a second companion. Nevertheless, it has 
the advantage of being simple, and it is worth noticing 
that all the SB in our sample satisfy this condition. 
Moreover, the planets that were not followed by RV ob- 
servations during a complete period are also discarded. 
Since this last condition results in removing the major- 
ity of the planets with periods longer than 2200 days, 
this value is adopted as a selection limit of the sample. 

— It seems that the planets found in binary or mul- 
tiple stellar systems have eccentricities smaller than 
the others when their periods are less than 30 days 
(Eggenberger et al. I2004J) . Three more planets are re- 
moved from the sample for this reason. 

— Multiple planetary systems are supposed to be differ- 
ent from those with single planets, since the eccentric- 
ities may be affected by resonant perturbations. It is 
worth noticing that the distribution of these systems 
around the median eccentricity is not significantly dif- 
ferent from that of the single planets : we find an excess 
of only one planet with a low eccentricity (when an ex- 
cess of large e is expected), and the probability to get 
this excess or a larger one just by chance is as large 
as 38 %. Nevertheless, for security, we still discard 14 
planets belonging to multiple planetary systems. 

A sample of 72 exoplanets orbiting main-sequence stars 
with periods shorter than 2200 days and with reliable or- 
bits is thus finally obtained. 

4.2. The exoplanets in the (P SI - e) diagram 

As for the SB, the planets are plotted on the (P sr - e) 
diagram (Fig. 0J) , in order to investigate the effects of tidal 
circularization. The sample is split into two nearly equal 
groups, one containing the planets with minimum mass 
less than 2 Jupiter, and one with the planets heavier than 
this limit. In contrast to the SB, for which circular orbits 
and moderate eccentricities are mixed in a small range of 
P sr , the separation between the circularized orbits and the 
others is remarkably well determined, at P cu toff = 5 days. 
The most eccentric planetary orbit, HD 80606b (Naef et al. 
2001}, is actually found for this period, but Wu & Murray 
(2003) demonstrated that it may be excited by a distant 
companion through the Kozai mechanism. 

It appears in Fig. 0] that the clustering of planets with 
P < 10 days (Udry et al. 120031) is even more marked when 
P sr is used in place of P : the planets are concentrated 
in orbits with the semilatus recta corresponding to P sr 
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Fig. 4. Same as Fig.|2lfor the planets. The circles are the 
planets with minimum mass < 2 Jupiter, and the open 
squares refer to the others. The thin dotted line is the 
limit esup corresponding to a maximum period of 2200 
days. 



0.8 



0.6 



0.4 



0.2 - 



0.0 



• •• 



10 



100 



1000 



Fig. 5. Same as Fig.|3for the exoplanets. The symbols are 
the same as in Fig.0] The full line is the median eccentric- 
ity of the planets with minimum mass < 2 Jupiter, and 
the dashes refer to the masses heavier than 2 Jupiter. The 
thin dotted line is the eccentricity eMax corresponding to 
a circularized period of 5 days. 



between 2.5 and 10 days, since we count 18 planets in this 
range, but only 1 between 10 and 20 days. 

It would be relevant to check if the eccentricities are 
correlated with P sr , but our sample does not permit this : 
the detection of the planets is far from complete, and the 
incompleteness increases with the period. Therefore, since 
the eccentric orbits correspond to a longer period for a 
fixed P sr , the (P sr - e) diagram of the planets is biased 
against large eccentricities. For that reason, the (P sr - e) 
diagram of the planets can be used only for investigating 
the circularization limit. As a consequence, it cannot be 
used to compare the SB to the exoplanets in Sect. El 

4.3. The exoplanets in the (P - e) diagram 

The (P - e) diagram of the exoplanets is given in Fig. [S] 
Only 1 planet above 2 Jupiter masses has P between 5 and 
70 days. This paucity of heavy-mass planets with short 
periods has already been pointed out (Zucker & Mazeh 
2002} Udry et al. 20O- and it makes the median eccen- 
tricity of these planets unreliable in this range of period. 
A two-sided test based on the common median for P in 
the range 5 to 2200 days shows that the probability to get 
differences at least as large as that obtained is 60 %. It 
is then quite possible that the eccentricities of planetary 
orbits are not related to the masses of the planets. This 
question is considered again in Sect. 

The median eccentricities of the exoplanets in Fig. [3] 
seem approximately constant. A test confirms this impres- 
sion : the Spearman coefficient of the 52 planets with 
P > 20 days and e < 0.78 is 0.21, providing a thresh- 
old between 10 and 20 %. It is thus not possible to rule 
out the hypothesis that, apart from circularization due to 
tidal effects, the distribution of the eccentricities of the 
planets is the same for any period between 20 and 2200 
days. Therefore, if the eccentricities are modified by mi- 



gration, they are changed almost independently of the pe- 
riods. However, this applies essentially to periods longer 
than 200 days, since we have very few planets between 20 
and 200 days. 

4.4. (P - e) diagram and metallicities 

Santos et al. (2004) have shown that the percentage of 
stars harboring planets jumps from less than 5 % to more 
than 20 % when stars with [Fe/H] larger than 0.2 are 
considered. In order to see if this limit of 0.2 dex also cor- 
responds to other orbital properties, we use the [Fe/H] of 
Santos et al. to distinguish the planets orbiting "metallic" 
stars and the others. In the (P - e) diagram, we count 27 
"metallic" planets, of which 14 have eccentricities smaller 
than the median. We conclude with a significance of 79 % 
that metallicity is not related to the eccentricity, confirm- 
ing the result obtained by Santos et al. (2003) with an- 
other test. 

5. The planets compared to the binaries 

5.1. The period-eccentricity diagram 

Non-twin binaries and planets are plotted in the (P - e) 
diagram in Fig. H3 The median eccentricities are derived 
from the systems with e < eMax- In order to make the 
comparison free of differences in the tidal circularization, 
eMax is derived from Eq. © assuming P cu toff = 10 days for 
both samples. It appears clearly that, although all objects 
are distributed in the same area of the (P - e) diagram, 
the planets have eccentricities that are on average smaller 
than those of the SB. The test of the distribution around 
the common median confirms this discrepancy : Among 53 
planets included in a sample of 102 objects, we count 33 
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Fig. 6. The exoplanets compared to the non-twin binaries 
in the (P - e) diagram. The circles refer to the planets, 
and the open squares to the SB with q < 0.8. The thin 
dotted line is eMax, the maximum eccentricity when the 
circularization limit is P C utoff = 10 days. The median ec- 
centricities are derived after discarding the systems above 
this limit, in order to have the same selection effects for 
both samples. The median eccentricity of the planets is 
given by a full line, and that of the SB by a a dashed line. 



planets below the median eccentricity. The null hypothesis 
is rejected at the 1.7 % level of significance. 

It appears from Fig. El that several planets with long 
periods and minimum masses below 2 Jupiter have very 
small eccentricities. Therefore, although we have seen that 
the (P - e) relation may be the same for all planets, it 
is relevant to compare the SB only to the planets with 
masses larger than 2 Jupiter. When the planets above 2 
Jupiter masses are compared with all the non-twin SB, 
the null hypothesis is rejected again at the 4.8 % level of 
significance, confirming that the heavy-mass planets have 
less eccentric orbits than the binaries. 

5.2. The intrinsic distribution of eccentricities 

A direct comparison of the distributions of the eccentric- 
ities of planets and of SB is not feasible, since the possi- 
ble range of eccentricities varies with the period, and the 
distribution of P is not the same for planets as for SB. 
Fortunately, another approach may be used to visually 
compare these objects, which is to derive the original dis- 
tributions of eccentricities corrected for the bias coming 
from tidal circularization. The low significance values of 
the Spearman correlation tests performed above allows us 
to assume that, apart from the area affected by circular- 
ization, the eccentricity distribution does not vary with 
the period. Therefore, it is possible to derive the intrin- 
sic distribution of the eccentricities for planets and for SB, 
using the method of the "nested boxes". For that purpose, 
we use Eq. (J2J to compute, for each system, the maximum 
eccentricity unaffected by tides, eMax- We then apply the 
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Fig. 7. Distribution of the eccentricities of SB and planets 
corrected for tidal circularization. 



method given in the Appendix. The results are shown in 
Pig. El 

The frequency of planets with low eccentricities is 
mainly due to the planets with masses below 2 Jupiter. 
The largest difference between the two groups of planets 
is obtained for e = 0.20 exactly : 40 % of the planets below 
2 Jupiter masses have orbits with e < 0.20, instead of 18 % 
for the others. However, a Smirnov test indicates that this 
discrepancy is far from sufficient to reject the hypothesis 
that all planets obey the same distribution, since the rejec- 
tion threshold is as large as 38 %. An anonymous referee 
wondered whether this test could be affected by a bias re- 
lated to the radial velocity semi-amplitude. However, this 
bias is unfavorable to the planets with the lowest masses, 
since they are the most difficult to detect, and its effi- 
ciency is maximum when they have eccentric orbits (see 
Sect. 15. t hereafter). As a consequence, correcting this bias 
would slightly decrease the proportion of e < 0.2 among 
the planets below 2 Jupiter, and then still increase the 
threshold of the test. We conclude again that we can as- 
sume that the eccentricity distribution of the planets does 
not depend on mass. 

The differences pointed out by the tests and by the 
comparison of the medians are clearly visible in this fig- 
ure. The maximum of the distribution is around 0.4 or 
0.5 for the non-twin SB, and between 0.3 and 0.4 for the 
planets. Moreover, the distribution of e decreases rapidly 
after the maximum for the planets, but it is rather flat 
until 0.8 for the non-twins SB. A comparison between the 
median e of the exoplanets and of the twins could suggest 
that these two kinds of objects have similar distributions 
of eccentricities. However, differences appear in Fig. : 
instead of exhibiting a maximum like the exoplanets, the 
distribution of eccentricities of twins is nearly flat over a 
wide range, from e = 0.1 to e = 0.7. 

5.3. Is this difference real ? 

We now want to check that the lack of planets with large 
eccentricities is not due to a selection effect against the 
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detection of these systems. In contrast to the SB, the de- 
tection of planets is far from complete, and our sample 
is obviously biased in favor of those which are easiest to 
detect. The most obvious bias is against the detection of 
long period systems, but this does not affect the reliabil- 
ity of our test based on the median e; it just decreases 
the contribution of the long-period planets. However, an- 
other bias is directly related to the eccentricity : a large 
eccentricity increases the semi-amplitude in RV, but, at 
the same time, it decreases the rms of the RV measure- 
ments. Therefore, the detection of a system close to the 
limit of the instrument is more difficult when the eccen- 
tricity is large. Another consequence of this effect is a bias 
in the distribution of the periastron longitude, to. The or- 
bits with u> around or 7r are more difficult to detect 
than those with u around tt/2 or 37r/2. This is visible, al- 
though not very significant, in our sample of planets: we 
count only 19 orbits with uj E [-7r/4, +7r/4] U [37r/4, 57r/4] 
among 43 planets with eMax > e > 0.1 (the orbits with 
e < 0.1 are not taken into account since u> is then not 
reliable). 

Simulations have been performed to investigate if this 
bias may explain the discrepancy between planets and SB. 
Each planet receives the eccentricity of a SB, randomly 
taken among the 10 SB with periods closest to that of the 
planet. The periastron longitude of the planet is randomly 
generated, and 3 radial velocity measurements are pro- 
duced for 3 epochs randomly chosen, adding errors drawn 
from the residual rms of the true orbit; (in reality, each 
star observed for planet detection receives much more than 
3 observations, but our aim is to derive an upper limit 
to the bias). When the standard deviation of the simu- 
lated RV is larger than the threshold corresponding to 
P(x 2 ) = 1 %, the planet is counted as detected; if not, 
another eccentricity is generated, and the simulation of 
the detection is performed again, until the planet satisfies 
the detection condition. When the complete sample has 
thus been detected by the simulation, the test of the me- 
dian eccentricity in the (P - e) diagram is performed. The 
simulation of the diagram is repeated 50 000 times. 

It appears from this calculation that the effect of the 
bias is to shift on average 0.6 more planets below the 
common median. Assuming that the number of planets 
below the median would be 32 in the absence of bias 
(instead of 33, see Sect. 15. If) , the rejection threshold of 
Hq becomes 4.7 %. This is still small enough to maintain 
rejection. At the same time, we count the planets with 
oj G [— 7r/4, +7r/4] U [37r/4, 57r/4] which are detected in the 
simulation. Among the planets with e > 0.1, their propor- 
tion is 45.5 %, in very nice agreement with the observed 
one, which is 19/43 = 44 %. We conclude then that the 
bias against detection of orbits with large eccentricities 
cannot explain the excess of planets with e smaller than 
the common median in the (P - e) diagram. 



6. Discussion and conclusion 

We have found some relevant features in comparing the 
eccentricities of the SB to those of the exoplanets : 

— The (P sr - e) diagram, based on P ST defined in Eq. Q , 
is a powerful tool for determining the limit between 
the circularized orbits and the others. The contrast 
between the two areas in the diagram suggests that 
the tidal effects are efficient only when the semilatus 
rectum of the orbit is less than a fixed limit. 

— The transition from the circularization to the part not 
affected by tidal effects looks sharper for the exoplan- 
ets than for the SB. For the SB, it corresponds to 
P sr between around 5 and 10 days, in agreement with 
Mathieu & Mazeh l|1988f> . Duquennoy et al. I|1992|l . 
Mathieu et al. I[1992fl and Witte & Savonije (2002), 
who consider that circularization is not restricted to 
the time of binary formation, but is still progressing 
during the whole lifetime of the main sequence compo- 
nents. For the exoplanets, the fast transition observed 
for planets with different ages is consistent with the 
idea that the tides were efficient only during the for- 
mation of the system. The inefficiency of tides for a 
formed planet is in agreement with the circularization 
time derived by Zahn l|1977[) . which is a function of 
(1 + q)/q. If the planets were brought closer to their 
host stars by migration, that means that migration oc- 
curred when the formation of the planets or of the host 
stars was not completed. Tidal circularization was then 
dominated by the tidal bulge on the planet, which was 
hotter, and therefore larger than it is today. 

— Before tidal effects had modified them, the eccentric- 
ities of planets or binaries were not strongly related 
to the periods, or to the angular momenta. It is even 
quite possible that they were not correlated at all with 
these parameters, since the absence of correlation is 
not clearly rejected by statistical tests, and also be- 
cause our assumption that tides did not affect at all 
the orbits with P sr > P cu toff is possibly a bit too sim- 
ple. Since the planets in the sample are supposed to 
have migrated, this suggests that migration did not 
alter the eccentricities significantly; alternatively, it is 
possible that the eccentricities were modified, but al- 
most independently of the final periods. 

— The exoplanets have orbits with eccentricities signifi- 
cantly smaller than those of the SB with the same pe- 
riod and with mass ratios larger than 0.8 (the non-twin 
binaries). A similar feature has already been pointed 
out for P < 50 days (Udrv l2T)01f) . but neglecting the 
difference between the distributions of periods of bi- 
naries and of planets. Moreover, it is now certain that 
the low eccentricities of planets are not an effect of 
the selection of the observed sample. Additionally, it 
seems that the distributions of the eccentricities are 
not related to the masses of the companions, neither 
among the non-twin binaries, nor among the planets. 
Therefore, this discrepancy is probably not an effect 
of the low masses of planets in a formation/evolution 
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process common to planets and binaries : this would 
imply a process depending on the secondary mass, but 
only around the transition between stellar and plane- 
tary companions. In fact, the SB most similar to the 
planets are the twins, perhaps because these systems 
were also interacting with a disk at the time of their 
formation, as Tokovinin (2000) suggested. 

Our most relevant conclusion is that the eccentricities 
of the exoplanet orbits are rather in favor of the hypoth- 
esis that exoplanets and binary stars are not the prod- 
ucts of the same physical process. After the "brown dwarf 
desert" (Halbwachs et ;i 1.127)00; separating the stellar com- 
ponents from the planets in the distribution of the sec- 
ondary masses, it is another argument in that sense which 
was derived from statistical investigations. 
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Appendix A: The method of the nested boxes 
applied to the distribution of eccentricities 

This method was initially dedicated to the derivation of 
a bias-free distribution of mass ratios of visual binaries 
(Halbwachs 1983 , Halbwachs et al. H997[l . as is explained in 
detail in Halbwachs (2001). It is adapted hereafter to the 
derivation of the intrinsic distribution of the eccentricities. 

We consider a sample with periods P larger than 
-Pcutoff, the period corresponding to tidal circularization, 
as explained in Sect. 12.11 For each system, P and Pcutoff' 
are used to derive eMax, the maximum eccentricity of the 
orbits unaffected by tidal effects. The systems having ec- 
centricities e > eMax j if any, are discarded from the sam- 
ple. 

The principle of method is as follows : 

— We define a first "box" by setting a small minimum 
value of eMaxi called e\ hereafter. The eccentricities of 
the systems having eMax > &i are used to derive a first 
estimation of the intrinsic distribution of eccentricities, 
/i(e), which is defined from e = 0uptoe = ei. Since e± 
is small, this distribution is based on a large number 
of systems, and it is fairly reliable; unfortunately, it 
concerns only a small range of eccentricities. 

— A second box is defined, using a limit ei > e%. The 
eccentricities of the systems with eMax > e-2 provide 
a second estimation of the distribution, /2(e), which 
applies to e < e2- We have then an estimation of /(e) 
which is valid between ei and e2, but for e < ei, /2(e) 
is less reliable than /i(e), since it is derived from fewer 
systems in this range. 



— The two boxes defined above are "nested", since the 
systems with e < ei belonging to the second box are 
all also present in the first box. This common part 
is used to derive /(e) by connecting / 2 (e > ei) to 
A(e < el). Let N\ and N 2 be the numbers of systems 
in the first box and in the second box, respectively. If 
the second box contains ri2 (ei) systems with e < e±, 
the best estimation of /(e) is : 

/1.2(e) = /i(e) + 6(e - ei )_^/ a ( e ) (A.l) 

«2(ei) 

where 6 is the Heaviside function; it is unnecessary to 
normalized the distributions so early in the calculation. 

— We may still add several boxes, using limits e^ with in- 
creasing values. If Ni-i is the norm of the distribution 
derived from the boxes 1 to i — 1, fi i—x, and if box 
i contains rij(ej_i) systems with e < e,_x, Eq. IjA.lfl 
becomes : 

/i,i(e) - / M -i(e) + 9(e - e^) /,(e) (A.2) 

ni(ei-i) 

The distribution is normalized after adding the last 
box. 

In practice, the values of the e^ terms are i x 0.01 when 
the eccentricities are provided with two decimals, in order 
to take into account all the systems with e < eMax- At the 
end of the calculation, the 0.01-bins are merged into 0.1- 
bins, in order to make the final distribution more readable. 
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